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Abstract. We report on the XMM-Newton observation of the Be/neutron star X-ray binary system RX J0146. 9+6121, a long 
period (~ 23 m) pulsar in the NGC 663 open cluster The X-ray luminosity decreased by a factor two compared to the last 
observation carried out in 1998, reaching a level of ~ 1 x 10'^'' erg s~^, the lowest ever observed in this source. The spectral 
analysis reveals the presence of a significant excess at low energies over the main power-law spectral component. The soft 
excess can be described by a black-body spectrum with a temperature of about 1 keV and an emitting region with a radius of ~ 
140 m. Although the current data do not permit to ascertain whether the soft excess is pulsed or not, its properties are consistent 
with emission from the neutron star polar cap. This is the third detection of a soft excess in a low luminosity (~ 1 x 10''^ erg 
s"'^) pulsar, the other being X Per and 3A 0535+262, suggesting that such spectral component, observed up to date in higher 
luminosity systems, is a rather common feature of accreting X-ray pulsars. The results on these three sources indicate that, in 
low luminosity systems, the soft excess tends to have a higher temperature and a smaller surface area than in the high luminosity 
ones. 
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1. Introduction 

High Mass X-Ray Binaries (HMXRBs) are binary systems 
consisting of a neutron star (NS) or, less frequently, a black 
hole, accreting matter from a high mass early type star. Based 
on the nature of the mass donor star, either a supergiant of O/B 
spectral type or a Be star, HMXRBs can be divided in two sub- 
groups showing different variability properties. HMXRBs with 
supergiant companions tend to be persistent sources, although 
some of them display high and low states with X-ray flux vary- 
ing by large factor (~ 100 or more). On the other hand, neu- 
tron star Be X-ray binaries (BeXRBs) are generally transient 
sources, owing to the long term variability of the equatorial 
discs suiTounding Be stars and/or the orbital eccentricity. For 
both subgroups, when the compact object is a NS, the X-ray 
spectra are well described by a rather flat power-law between 
0.1 and 10 keV (photon index ^^1) followed by a high-energy 
cutoff. 

With the advent of imaging X-ray satellites a large num- 
ber of BeXRB sy stems has been discovered in the Small 
MageUanic Cloud faaberl & Sasakilj2o'oo[ llsrael et al.lj200C ; 



lYokogawa et al.l2003t iMacomb et al fcool ISasaki et al 12003 ; 
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IHaberl & PietschI l2004l) . The observation of these sources, 
unaffected by the high interstellar absorption present in the 
Galactic plane, makes it possible to study in detail their X- 
ray spectra extending down to energies of a few hundred 
eV. This has allowed to discover that most of them have a 
marked soft excess abo ve the power-law model (■Nagase.2002i: 
iHaberl & Pietschl2005l) . 

RX J0146.9+6121 is a BeXRB hosting a neutron star char- 
acterized by a rotational period of about 23 minutes, among the 
longest observed in X-ray pulsars. The pulsations were discov- 
ered with non-imaging instruments on boar d EXOS AT and ini - 
tially attributed to a different nearby source dWhite et alJl987l). 
Subsequent observations clarified the picture (Meregh etti etalJ 
1993: Hellier 1994) and led to the optical identificat ion of RX 
J0146.9+6121 wifli flie BO Ille star LS I +61° 235 JCoe et alJ 
1993). This st ar i s a member of the open clus ter NGC 663 
(iTapJa et al.i iil99ll IPabregat et alJ il996: Pigulski et alJ 1200 ih 
for which a d istance in the range between 2 and 2.5 kpc has 
been derived ("Kharchen ko et al.ll2005t IPandev et al'] l2005g). In 
the following we adopt d = 2.5 kpc. Thanks to its relatively 
small distance, RX J0146. 9+6121 is not much absorbed and is 
therefore a good target to investigate the properties of the soft 
X-ray emission in a Galactic source. Here we present the re- 
sults of a recent observation obtained with the XMM-Newton 
satellite, providing the most sensitive observation of this source 
ever obtained below 2 keV. 
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2. Observations and data reduction 

RX J0146.9+6121 was observed by XMM-Newton between 
22:40 UT of 2004 January 14 and 10:20 UT of 2004 January 
15. Since the main target was the NGC 663 open cluster, RX 
J0146. 9+6121 was detected at an off-axis angle of 9.3 ^ All the 
three EPIC instruments, i.e. the pn ca mera llStriider et alJ200lh 
and the MOSl and M0S2 cameras dTurneret alJl200ll) . were 
active and operated in Full Frame mode. For all of them, the 
Medium thickness filter was used. 

We used version 6. 1 of the XMM-Newton Science Analysis 
System (SAS) to process the event files. After the standard 
pipeline processing of the data, we looked for possible peri- 
ods of high instrumental background, due to flares of soft pro- 
tons with energies less than a few hundred keV. We found that 
the first ~ 15 ks of the observation were affected by a soft- 
proton contamination. However, since RX J0146. 9+6121 was 
detected with a count-rate (~ 1 cts s^^) much higher than that 
of the background, the soft-protons during the bad time inter- 
vals have a negligible effect on the source spectral and timing 
analysis (their count rate in the source extraction area is less 
than 0.01 cts s^^). Hence we used for our analysis the data of 
the whole observation, corresponding to exposure times of 35.4 
and 41.2 ks in the pn and MOS cameras, respectively. 

3. Timing anaiysis 

During our observation some flux changes on ~ hour timescale 
were present, as shown by the background subtracted light 
curve plotted in Fig. [2 This is based on the data in the 0. 15-10 
keV range obtained from the three EPIC cameras. A bin size 
of 1 .4 ks, corresponding to one spin period, has been chosen to 
avoid the effects due to the periodic pulsations. Variations up to 
~ 20 % around the average level of ~ 1.9 cts s^^ are evident. 
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Fig. 1. Background subtracted light 
J0146.9+6121 in the energy range 0.15- 
time bin of 1.4 ks (i.e. 1 pulse period). 



curve of RX 
10 keV, with a 



To obtain a measure of the pulse period, we converted the 
times of arrival to the solar system barycenter and performed a 
folding analysis using the source events of three cameras. By 
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Fig. 2. Background subtracted light curves of RX 
J0146.9+6121 in the energy ranges 0.3-2, 2-4.5, 4.5-7 
and 7-10 keV, folded at the best-fit period of 1396.14 s. The 
vertical lines indicate the phase intervals used for the spectral 
analysis. 

fitting the x^ versus trial period curve with the appropriate sine 
function as described in ' Leah vl ( Il987l) . we derived a period of 
1396.14 ± 0.25 s. In Fig. |2]we show the folded light curves 
in four energy intervals (0.3-2, 2^.5, 4.5-7 and 7-10 keV) 
chosen to allow a direct compariso n with the hght curves previ- 
ously derived with RossiXTE data (iMereghetti et al..2000.) . The 
pulse profile, characterized by a broad peak, is clearly energy- 
dependent, with the maximum shifting from phase ^ 0.4-0.5 
at E < 2 keV to phase - 0.6 at E > 7 keV. 

The hardness ratio (HR) between the light curves above and 
below 2 keV, reported in Fig.|3] is characterized by a slow and 
regular increase from its minimum to its maximum value, fol- 
lowed by a rather irregular decrease. This plot, where we have 
divided the pulse period in ten phase intervals, shows that there 
is not a simple correlation between the hardness and the total 
count rate: we observe the same HR value at completely dif- 
ferent count rate levels, but also very different HR values for 
the same count rate. This indicates that the spectral hardness of 
RX J0146.9+6121 does not depend in a simple way on its flux 
level. 

4. Spectrai anaiysis 

For the source spectra we used an extraction radius of 30" 
around the source position in the case of the M0S2 and pn cam- 
eras; for the MOSl camera the extraction radius was reduced 
to 15", since the source was imaged close to a CCD gap. We 
checked with the SAS task epatplot that no event pile-up af- 
fected our data, then we accumulated all the events with pat- 
tern range 0-4 (i.e. mono- and bi-pixel events) and 0-12 (i.e. 
from 1 to 4 pixel events) for the pn and the two MOS cam- 
eras, respectively. The background spectra were accumulated 
on large circular areas with no sources and radius of 210" and 
120" for the MOSl and the M0S2 camera, respectively. For the 
pn camera the background spectrum was extracted from a cir- 
cular region of the same radius and at the same CCD rows of 
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Fig. 3. Hardness ratio of the two light curves of RX 
J0146. 9+6121 in the energy ranges 0.3-2 and 2-10 keV, as a 
function of the pulse phase. The HR values are obtained by 
folding the light curves at the best-fit period of 1396.14 s. The 
pulse period is divided in 10 phase bins, which are identified 
by progressive numbers. 
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Fig. 4. Top panel: pn spectrum of RX J0146.9+6121 with the 
best-fit power-law {dashed line) plus black-body (dotted line) 
model. Middle panel: residuals (in units of a) between data 
and model in the case of the single power-law. Bottom panel: 
residuals in the case of the power-law plus black-body. 



the source region. We generated ad hoc response matrices and 
ancillary files using the SAS tasks rmfgen and arfgen. All spec- 
tra were rebinned with a minimum of 30 counts per bin and 
fitted in the energy range 0.3-10 keV using XSPEC 1 1 .3.2. 

After checking that separate fits of the three spectra gave 
consistent results, we fitted the spectra from the three cam- 
eras simultaneously in order to increase the count statistics and 
to reduce the uncertainties. The fit with an absorbed power- 
law yielded Nr = (6.24 ± 0.14) x lO^^ cm^^ j^^^j pjjQ. 
ton index F = 1.36 ± 0.02, but with large residuals and 
X^/d.o.f. = 1.283/1480'. The addition of a blackbody compo- 
nent improved the fit quality significantly (Fig.|4}: we obtained 



^21 



^ r = 1.34; 



-0.05 



and kTsB 



Nh = (5.09l°J4) X 10^ 
l.lltg;^^ keV, with x^/d.o.f. = 1.036/1478. The emission sur- 
face of the thermal component has a radius Rbb = 140^14 m 
(for d = 2.5 kpc). The unabsorbed flux in the energy range 0.3- 
lOkeVis/x ~ 2x 10~" ergcm^^ g-i ,^^0^24 % of which 
is due to the blackbody component. 

We also attempted to fit the soft part of the spectra with 
other emission models, such as mekal, thermal bremsstrahlung 
and broken power-law. In all these cases the results were worse 
than those obtained with the blackbody model, since the re- 
duced chi-squared was higher, the residuals were larger and/or 
the best-fit parameter values were unrealistic. 

We did not find significant evidence for emission lines. By 
adding Gaussian components at various energies to the above 
model, we estimated upper limits of the order of 150 eV for the 
equivalent width of lines in the 6-7 keV energy range. 

5. Phase-resolved spectroscopy 

Prompted by the results described above we performed a 
phase-resolved spectroscopy in order to study in more detail 



the source behavior We defined phase by matching the folded 
light curve to that observed in 1996 with RossiXTE and ex- 
tracted the background subtracted spectra for the same four 
phase intervals used in Mereahetti et al. (2000) and indicated 
in Fig|2] The first step was to fit all of them with the best- 
fit power law plus blackbody model of the phase averaged 
spectrum, leaving only the relative normalization factors free 
to vary. The ratios of the four spectra of each instrument to 
these renormalized average models show significant residuals 
and clearly demonstrate the spectral variability as a function of 
the pulse phase. 

We then fitted the four spectra independently. In all cases, 
the absorbed power-law model was not satisfactory, while the 
addition of a black-body component significantly improved the 
fit. Therefore we used this model for all the spectra, leaving all 
the parameters free to vary; the results are reported in Tab.[l] 

Table 1. Best-fit spectral parameters for the phase-resolved 
spectroscopy of RX J0146. 9+6121, in the case of the indepen- 
dent fit of the four spectra. Errors are at 90 % confidence level 
for a single interesting parameter. Nh and kTsB are measured 
in units of lO^'^ cm^^ and keV, respectively. 



Spectral 




Phase Interval 




Parameter 


A 


B 


C 


D 


Nh 

r 

klBB 


4 8+"-* 

^■0-0.7 

, 4y+0.19 

J--'l'-0.20 


5 4+0.4 

0-*-0.3 

1 oq+0.09 

1 n'=;+°i° 


1 qi+0.35 

J^"JJ^-0.23 
1 34+0.17 
J^-J*-0.18 




/tot 


1.81 


3.22 


1.74 


1.00 


/pL 


1.34 


2.52 


1.10 


0.70 




(74 %) 


(78 %) 


(63 %) 


(70 %) 


/bb 


0.47 


0.70 


0.64 


0.30 




(26 %) 


(22 %) 


(37 %) 


(30 %) 


Xi/d.o.f. 


0.972/308 


1.021/992 


1.085/352 


1.131/328 



Errors are at 90 % confidence level for a single interesting param- " Unabsorbed flux in the energy range 0.3-10 keV, in units of 10" 



eter 



erg cm 
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In order to investigate the relative variations of the two 
components with the period phase, we also fitted simultane- 
ously the four spectra forcing common values for Nh, F and 
kTsB- In this case we obtained Nh = (5.8±0.1) x 10^^ cm^^, 
TpL = 1.60 ± 0.02 and kTsB = l-36+oo3 keV, with x^/d.o.f. 
= 1.117/1989; the corresponding normalization values are re- 
ported in Tab. |2] In this interpretation the spectral changes as 
a function of the phase are reproduced by the variations in the 
relative contribution of the two components. 



Table 2. Best-fit values for the black-body and power-law nor- 
malizations, when the four spectra are fitted simultaneously 
with common values of Nh ((5.8 ± 0.1) x 10^^ cm~^), T 
(1.60 ± 0.02) and kTeB (l.SGtaog keV). EiTors are at 90 % 
confidence level for a single interesting parameter 



Spectral 




Phase Interval 




Parameter 


A 


B 


C 


D 


pa 


76 ±3 


130 ±2 


121 ±2 


87 ±2 




2 24+°-"'' 


q 9t-+0.05 


1 09+0.03 


o.89+°:S 


/tot 


2.04 


3.42 


1.77 


1.06 


/pL 


1.69 


2.46 


1.00 


0.67 




(83 %) 


(72 %) 


(56 %) 


(63 %) 


/bb 


0.35 


0.96 


0.77 


0.39 




(17 %) 


(28 %) 


(44 %) 


(37 %) 



" Radius of the blackbody component (in metres) for a source distance 

of 2.5 kpc. 

' Intensity of the power-law component in units of 10^'^ ph cm^^ s^^ 

keV"^ at 1 keV 

■^ Unabsorbed flux in the energy range 0.3-10 keV, in units of 10"^^ 

erg cm~^ s~^ 



However, the above results do not necessarily imply that 
the soft component is pulsed. In fact an acceptable fit (x^/d.o.f. 
= 1 .097/1989) can also be obtained by imposing that the black- 
body component parameters be the same in the four spectra, as 
shown in Tab|3] 



The above results are summarized in Fig.|5l where we show 
the pulse-phase dependence of the black-body temperature, 
the power-law photon index and the unabsorbed flux of the 
two components for both fits. Since they have a similar statis- 
tical quaUty, we can neither confirm nor deny that the thermal 
component is variable. 
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Table 3. Best-fit values for the power-law parameters, when 
the four spectra are fitted simultaneously with common values 
of Nh ((6.1I0:?) X 10^^ cm-2), kTeB (1.40l^;[54 keV) and 
Rbb (IIII4 m). Errors are at 90 % confidence level for a single 
interesting parameter 



Fig. 5. Pulse-phase dependence of the black-body tempera- 
ture, the power-law photon index and the unabsorbed flux of 
the two components (in units of 10^^^ erg cm^-^ s^^), in the 
case of common temperature and index (top) and of common 
black-body (bottom) 
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Parameter 


A 


B 


C 
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r 

Ttt 


1 ga+0-U3 

-'-•y°-0.04 

2 40+0'^^ 


1 ry+U.Ol^ 
J^-O'-O.02 
47+0.05 


1 c:fi+"'"* 
J^.OD_0.03 
1 oq+0.03 


9 ^o+0.i3 
07+0-04 


/tot 


2.08 


3.43 


1.81 


1.19 


rb 
/pL 


1.36 


2.71 


1.09 


0.47 




(65 %) 


(79 %) 


(60 %) 


(40 %) 


Jbb 


0.72 


0.72 


0.72 


0.72 




(35 %) 


(21 %) 


(40 %) 


(60 %) 



" Intensity of the power-law component in units of 10 ^ ph cm ^ 
s-i keV"^ at 1 keV 

' Unabsorbed flux in the energy range 0.3-10 keV, in units of 10~^^ 
erg cm~^ s~^ 



6. Discussion 

The XMM-Newton data described here were obtained almost 6 
years after the latest X-ray observation of RX J0146.9H-6121, 
which was performed by BeppoSAX on 1998 February 3 
( Mereghetti et al. 2000) . Therefore it is interesting to compare 
our results with those obtained in the past. In Fig. |6|we show 
the long term evolution of the source spin period and lumi- 
nosity since the time of the EXOSAT discovery in 1984. If we 
exclude the first observation, when the source was in outburst, 
a linear fit of all the periods gives a spin-up at an average rate 
of P = — (4.6io'2) X 10^* ss^^, similar to that measured until 
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1998. This result suggests that during the 6 years between the 
BeppoSAX and the XMM-Newton observations, the momen- 
tum transfer onto the neutron star has proceeded with no major 
changes. The flux detected by XMM-Newton corresponds to a 
luminosity of ^ 1.5 x 10'^'* erg s^^, a factor ^ 2 lower than 
the minimum level observed in the p revious years , indicating 
that, after the outburst of July 1997 faaberlet alJll998.) . RX 
J0146. 9+6121 has been continuously fading. 

Comparison wi t h the RossiXTE results reported in Fig. 2 
of Mereghetti et alJ OOOOl) shows that, despite the luminosity 
variation, the shape of the pulse profiles in the energy interval 
common to the two instruments (2-10 keV) has not changed. 
The XMM-Newton data also confirm that, above 2 keV, there 
is a significant spectral softening in the initial rising part of the 
pulse (phase interval A). 



1400 




1985 1990 1995 2000 

Date (y) 

Fig. 6. Pulse period and luminosity history of RX 
J0146.9+6121 from 1984 August 27 to 2004 January 15. 
The solid line in the upper panel is the best-fit linear model 
of the reported values and has a slope P ~ —4.6 x 10^® 
s s~^. The luminosities of the lower panel are based on the 
unabsorbed flux in the 2-10 keV energy range and on a source 
distance of 2.5 kpc. 



The main advantage of the data reported here, with re- 
spect to previous observations of this source, is the better cov- 
erage of the energy range below 2 keV. In comparison with 
the BeppoSAX results, the phase averaged spectroscopy gives 



a smaller hydrogen column density: Nh 



(5.09 



+0.24N 
-0.23-' 



10^1 cm-2 instead of (1.2±0.3)xl022 cm^^ jjjis value is 
in better agreement with that expected fro m the o ptical ob- 
servations which give E{B - V) = 0.93 dReig et al. 1997) . 
Assuming Ay = 3.1 E{B — V) and the average relation Ay 
= Nh X 5.59 x lO^ffcm"^ between optical extinction and 



X-ray absorption (^Predehl & Schmitt 1995) . this would predict 
Nh = 5.16 X lO^i cm-2. 

Moreover, we find evidence of a thermal component which 
was never observed before in this source. In the last years, sim- 
ilar soft excesses have been detected in several high mass X- 
ray binary pulsars. The properties of these sources are sum- 
marized in Table |4] Hickox et al. (:2004i) showed that a soft 



excess has been detected in all the X-ray pulsars with a suf- 
ficiently high flux and small absorption. In fact, most of the 
soft excess sources are at small distances and/or away from the 
Galactic plane (most of them are in the Magellanic Clouds). 
This suggests that the presence of a soft spectral component 
could be a very common, if not an ubiquitous, feature intrinsic 
to X-ray pulsars. These authors concluded that the origin of the 
soft component is related to the source total luminosity. When 
ix > 10'^* erg s~^ the luminosity and the shape of the soft 
component can be explained only by reprocessing of hard X- 
rays from the neutron star by optically thick accreting material, 
most likely near the inner edge of the accretion disk. In less 
luminous sources, with Lx < 10"^® erg s^^, the soft excess can 
be due to other processes, such as emission by photo-ionized 
or collisionally heated diffuse gas or thermal emission from the 
surface of the neutron star Finally, in the sources of intermedi- 
ate luminosity, either or both of these types of emission can be 
present. 
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Fig. 7. Total X-ray luminosity of the sources of Tab. |3 as a 
function of the interstellar absorption. 

In Fig. we have plotted the X-ray luminosity versus in- 
terstellar absorption of the pulsars with a soft excess. We can 
see that RX J0146.9+6121, together with X Persei and 3A 
0535+262, is at least one order of magnitude less luminous 
than all the other pulsars displaying a soft excess. We note that 
the ratio between the unabsorbed fluxes of the thermal and the 
power-law components in RX J0146. 9+6121 agrees well with 
the average value measured in the other sources. However, un- 
like some of the other sources, the soft excess could be fitted 
only with a black-body component, while any other emission 
model was rejected. 



Table 4. Orbital and spectral parameters forXBPs with a detected soft excess (SE). 
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(kpc) 
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(ergs s^^, keV) 


Flux 
(ei-gscm^^ s"^) 


(10^1 cm-2 


Lse/Lx 
) 




SE model" 


kTBB 

(keV) 


SE Pulses 


e 




Her X-I ^ 


Galaxy 


~5 


A7V 1.24 


1.0 X 10-^' (0.3-10) 


3.3 X 10-'^ 


0.05 


0.04-0.10 




BB. BB-hLE 


0.09-0.12 
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SMCX-1^ 


SMC 


65 


BO lb 0.7 


2.4 X 10^** (0.7-10) 


4.7 X 10"" 


2-5 


0.036 




BB. TB, SPL 


0.15-0.18 
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LMC X-A^ 


LMC 


50 


07 UI-O IV 13.5 


1.2 X 10^** (0.7-10) 


4.0 X 10"" 


~0.5 
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0.15 
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1.9 


B0.5 lb 283 


2.2 X 10^*= (2-10) 


5.1 X 10-** 


4.2 


--^0.01 




TB 




No 




y, 


XPer» 


Galaxy 


0.95 


09.5pe 837 


1.8 X 10^" (0.3-10) 


1.7 X 10-" 


1.5 


0.24 




BB 


1.45 


? 




r 


EXO 053109-6609.2" 


' LMC 


50 


B0.7 Ve 13.7 


4.6 X 10^^ (0.2-10) 


1.5 X 10-" 


6.9 


7 




MEK-fPL 




Yes 




A 0538-66" 


LMC 


50 


B2 Ille 0.069 


4.0 X 10^^ (0.1-2.4) 


1.3 X IQ-" 


0.8 


7 




BB,TB 


~0.2 


? 




1 


RXJOO47.3-731212 


SMC 


65 


B2e 263 


1.5,2 X 10^'' (0.7-10) 


3.0,4.0 X 10-^^ 


0.96,3.4 


0.03-0.09,0.68 




BB 


0.6, 2.2 


Yes 




RXJ0101.3-721li^ 


SMC 


65 


Be 452 


1.6 X 10^** (0.3-10) 


3.2 X IQ-" 


0.6 


7 




MEK 




7 




RXJ0103.6-72011'' 


SMC 


65 


05 Ve 1323 


0.8 - 7.5 X 10^'^ (0.2-10) 


1.6- 14.8 X 10-^^ 


1.9 


7 




MEK 




No 




s 


AX J0049.5-7323^^ 


SMC 


65 


B2Ve 751 


7 - 9 X 10^^ (0.2-10) 


1.4- 1.8 X IQ-^^ 


3.6 


7 




7 




Yes 




Rp 


AX J0058-720^^ 


SMC 


65 


Be 281 


1.2 X 10^^ (0.3-10) 


2.4 X IQ-" 


0.6 


7 




7 




Yes 




AXJ0103-722^^ 


SMC 


65 


BO IV 342 


1.8 X lO^'^ (0.3-10) 


3.6 X 10-" 


0.6 


7 




MEK 




? 




k 

at 


3A 0535+2621'^ 


Galaxy 


2.0 


09.7 Ille 103.4 


3.9 X 10^^ (2-10) 


8.2 X 10-^^ 


6.0 


0.35 




BB 


1.33 


? 




RXJ0146.9+6121 


Galaxy 


I^^ 


BOnie 1395 


1.5 X 10^" (0.3-10) 


2.0 X 10-" 


5.1 


0.25 




BB 


1.11 


? 




3 


" References: (l)'d£ 


il Fiume et al." 

(2000b); (5)1 
lJl200lh:tlO^ 


(■1998'), 


'Endo et al.'('2000);Ramsav et al.'('2002"); (2)'W( 


00 et al.'('1995);Paul et al.'('2002'): 
et al. (2000); (8) Haberl ( 1994), ( 
3); (12) The two set of values are 


;(3)'Wooetal.'('l996). 
Drlandinietal. (1998), 


,'La Barbera et al. 


'('200l'),'Naik&Paul 


1^2004'); 




(4) Yokoeawa et al. 


S;ohno et al. (2000); (6) Schulz et al. (2001); (7) Burderi 


Ki-evkenbohm et al. (2002); 


(9) di Salvo et al. 
respectively; (13) 




il998h.lCoburnetal 


Haberl et al. ( 2003); (11) Mavromatakis & Haberl ( 199 
al..k2003,).,Haberl & Pietscln2005,): (15),Yokosawa et ; 


base on Ueno et al. (2004) and Maiid^ 


:alJJ2004h. 


1 


ISasakietalJj2003h: 


tl4lSasaki et 


alJj2000al).lHaberl & PietsclJi200 


i): (16),Mukheriee & Paul,k2005,). 


*" For each source also the reference energy range is reported. 




















^ 


'^ The reported values refer only to the interstellar absorption, not to the source intrinsic absorption. 
















? 


'^ Spectral models used for the soft excess are: BB = blackbody: 


; TB = thermal bremsstrahlung; SPL = soft 


power-law or broken power- 


-law; MEK = MEKAL thin thermal model; 


COM = 


s 


Comptonization model; LE = broad low-energy line emission. Commas indicate separate 


fits, plus signs indicate fits with two components. 












8- 


" Pulsation of the emission component that traces the soft excess. 




















n. 



+ 
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Based on the results obtained bv lHickox et alJ yOOJ), the 
luminosity of RX J0146. 9+6121 is too small for an interpre- 
tation of the soft excess in terms of reprocessing of the hard 
X-ray photons in optically thick accreting material. Moreover, 
both the thermal emission and the hard X-ray reprocessing in a 
diffuse, optically thin gas around the neutron star are unlikely, 
since these processes would not give a black-body spectrum. 
We therefore favor the interpretation of the soft excess in RX 
J0146. 9+6121 as thermal emission from the neutron star po- 
lar cap. If we assume that the source is in the 'accretor' status, 
with matter accretion on the NS surface, the blackbody emit- 
ting radius of ~ 140 m is consistent with the expected size of 
the polar cap. In fact, if Mns = 1-4 Mq and Rns = 10^ cm, the 
source luminosity of ~ 10^^ erg s^^ implies an accretion rate 
M '^ 5 X 10^^ g s^^ and, adopting Bn s = 10^^ G, a m agne- 
tospheric radius R^ ~ 2.4 x 10^ cm JCampana et alJ[l998.) . 
In this case, based on the relation Rcoz ~ Rns (RNs/Rm)*^'^ 
JHickox et al. '2004'). we would obtain Rcoi ^ 200 m. If this de- 
scription is correct, we would expect to observe some variabil- 
ity of the soft component along the pulse phase. From this point 
of view the phase-resolved spectroscopy provides no conclu- 
sive results, since it proves that both a variable and a constant 
thermal component can account for the observed spectral vari- 
ability. 

We finally note that RX J0146.9+6121 is in many respects 
very similar to X Per and 3A 0535+262, which are very low 
luminosity Be/NS binary with a long pulse period. The soft ex- 
cesses of these sources have similar properties, since they tend 
to have a higher temperature (> 1 keV) and a smaller emission 
radius (^ 0.1 km) compared to the soft excesses observed in 
high luminosity systems, which have a temperature of about 0. 1 
ke V and an emission radius of a few hundred km. In these three 
low luminosity systems the soft excess contributes for 25-35 % 
of the total flux cCoburn et al...200L Mukheriee & PauL2005.) . 
Also in the case of X Per and 3A 0535+262 this excess has 
been attributed to the emission from the polar caps. 

7. Conclusions 

We have reported on the analysis of the data collected by 
XMM-Newton in a ^ 42 ks observation of the Be/neutron-star 
X-ray pulsar RX J0146.9+6121. 

The unabsorbed flux corresponds to a source luminosity 
Lx ~ 1 X 10'^"' erg s^^ in the 2-10 keV energy range, about 
50 % smaller than the lowest level ever observed from this 
source, indicating a monotonic source fading over long time 
scales. 

Thanks to the high effective area of XMM-Newton also 
at low energies, we could perform the first accurate spectral 
study below 2 keV for this source. In the phase-averaged spec- 
trum we have revealed the presence of a significant soft excess 
over the primary power-law component: this excess can be de- 
scribed by a black-body with kTsB ~ 1 keV, while any attempt 
to fit it with a different emission model was unsuccessful. 

The phase-resolved spectroscopy has confirmed the large 
spectral variability along the pulse period already observed 
above 2 keV. Unfortunately, with the current data, it is not pos- 
sible to derive compelling results on the phase variability of the 



soft excess component. Although the emission below 2 keV is 
clearly pulsed and the low energy part of the spectrum varies 
with the phase, we have shown that such variations can be ex- 
plained equally well by changes in the blackbody component 
or in the power law component alone. 

Clearly the relatively small distance and low interstellar ab- 
sorption toward RX J0146.9+6121 played a role in the possi- 
bility of detecting a soft excess in such a low luminosity pul- 
sar Comparison with other X-ray binary pulsars shows that, 
so far, the soft excess had been detected only in much brighter 
sources. The data reported here support the hypothesis that a 
soft thermal component is an ubiquitous emission feature of 
this class of sources. In this sense, it would be very interesting 
to use the large collecting area of XMM-Newton in long ob- 
servations of the faintest and longest period Be binaries, both 
in the SMC and in the Milky Way, such as, for example, the 
persistent low luminosity systems R X J0440. 9+443 1 and RX 
J1037.5-5647 JReig & Rochell999t) . 
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